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Two-dimensional (2D) topological crystalline insulators (TCIs) were recently predicted in thin 
films of the SnTe class of IV-VI semiconductors, which can host metallic edge states protected by 
mirror symmetry. As thickness decreases, quantum confinement effect will increase and surpass 
the inverted gap below a critical thickness, turning TCIs into normal insulators. Surprisingly, 
based on first-principles calculations, here we demonstrate that (001) monolayers of rocksalt IV-VI 
semiconductors XY (X=Ge, Sn, Pb and Y= S, Se, Te) are 2D TCIs with the fundamental band gap 
as large as 260 meV in monolayer PbTe, providing a materials platform for realizing two-dimensional 
Dirac fermion systems with tunable band gap. This unexpected nontrivial topological phase stems 
from the strong crystal field effect in the monolayer, which lifts the degeneracy between p x , y and 
p z orbitals and leads to band inversion between cation p z and anion p XiV orbitals. This crystal field 
effect induced topological phase offers a new strategy to find and design other atomically thin 2D 
topological materials. 


Topological insulators (TIs) have attracted extensive 
attention owing to their fundamental theoretical interests 
and potential applications m ■ The interplay between 
crystallography and electronic band topology has further 
given birth to a new type of topological phases, termed 
topological crystalline insulators (TCIs) [J[[5J- The first 
example of TCIs was recently predicted in SnTe class 
of IV-VI semiconductors [6j. They host even number 
of Dirac cones on their (001) and (111) surface [MS, 
which were subsequently observed mm- The metallic 
boundary states of TCIs are protected by crystal symme¬ 
try rather than time-reversal in TIs, and therefore these 
states can acquire a band gap under perturbations that 
break the crystal symmetry by either distortion [B1 ITB1 - H5] 
or electric field m ■ In addition to fundamental interests, 
TCIs may realize novel device applications 12D1 - !2!H . Re¬ 
cent studies have further predicted/proposed a wide class 
of new TCI materials in antiperovskites [24], pyrochlore 
iridates [25] , multi-layer graphene [26] and heavy-fermion 
compounds mm- 

More recently, based on tight-binding (TB) calcu¬ 
lations [25] . we predicted that two-dimensional (2D) 
TCIs exist in (001) multi-layer thin films of SnTe and 
Pbi_ a ,Sn a; Se(Te) [19] . It was found that thin films thicker 
than a critical thickness are topological nontrivial with 
metallic edge states protected by mirror symmetry z —>■ 
—z (z is normal to the film). As thickness of TCI thin 
films decreases, the quantum confinement effect will in¬ 
crease and eventually overcome the inverted gap, turning 
these 2D TCIs into normal insulators [T5]. On the other 
hand, it is known that when materials become atomically 
thin, the reduced dimensionality may lead to physical 
properties dramatically different from their 3D counter¬ 
parts, as found in graphene [So] EU, monolayer M 0 S 2 
[32l:3H], monolayer a-Sn [34] and FeSe !35J. At the same 
time, there are great interests in the exploration of nanos¬ 


tructure of TCIs in experiments [151 fHHMT] . These moti¬ 
vates us to study IV-VI semiconductors in the monolayer 
form, whose electronic structures may not be captured by 
the bulk TB parameters [25] used in previous works. 

Here, using first-principles calculations, we systemati¬ 
cally studied the electronic structures of (001) monolayer 
IV-VI semiconductors XY (X= Ge, Sn, Pb and Y= S, Se, 
Te) in the rocksalt structure. Remarkably, we found all 
of them are 2D TCIs and possess metallic edge states 
protected by z —t — z mirror symmetry, contrary to the 
expectation from the quantum confinement effect. This 
unexpected nontrivial topological phase originates from 
the band inversion due to the strong crystal field effect 
on p orbitals in the monolayer, which significantly low¬ 
ers (raises) the onsite energy of cation (anion) p z orbital, 
thereby facilitating or even driving the band inversion to¬ 
wards the TCI phase. It is worth noting that this novel 
mechanism does not rely on the electronic topology of the 
3D bulk states: for example, despite that PbTe is topo¬ 
logically trivial in the bulk form, monolayer PbTe is a 2D 
TCI with a large fundamental gap of 260 meV. Therefore 
our work reveals a unique advantage of atomically thin 
2D materials in realizing 2D topological phases. 


RESULTS AND DISCUSSION 

The band structure plot displayed in Fig. B a ) shows 
that monolayer SnTe is a semiconductor with a small in¬ 
direct gap located around the X point. Moreover, the 
nearby IV-shape conduction band and M-shape valence 
band provide a clear sign of band inversion at X point 
with a large inverted gap (denoted by E*) of about 
—0.6 eV at its optimized equilibrium lattice constant of 
L op t = 6.2 A. Negative sign of indicates an inverted 
gap. To confirm this band inversion, we performed first- 
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FIG. 1. (a) Band structure of monolayer SnTe, where the size 
of red dots indicates the weight of Sn p z orbital; the inset is 
the 2D Billouin zone of (001) thin films; (b) Band gap at X 
point (Ef) and fundamental gap ( E g ) of monolayer SnTe as 
a function of lattice constant ( L ). As L increases, monolayer 
SnTe undergoes a topological phase transition from \Nm\ = 2 
to Nm = 0. When L > 6.2 A, it is a fully gapped semicon¬ 
ductor, otherwise, it is a semi-metal; The inset is the band 
structure evolution around the transition point (red:L=6.4 A; 
blue: L= 6.6 A; black: L= 6.8 A). 


principles calculations of monolayer SnTe under a series 
of biaxial elastic strain mmm- As shown in Fig. 0b), 
when lattice constant ( L ) increases, E* monotonically 
enlarges and reaches zero at a critical lattice constant of 
L c = 6.6 A. Beyond L c , E* becomes positive, indicat¬ 
ing a trivial topological phase. The monotonical change 
of E* is very similar to those in Bi2Se3-type topolog¬ 
ical insulator [42], which comes from the different par¬ 
ities of bottom conduction band and top valence band 
at X points: one corresponds to bonding states and the 
other corresponds to anti-bonding states, whose energy 
changes oppositely under the biaxial strain. However, 



FIG. 2. Edge states of monolayer SnTe along [001] direction 
with lattice constant L = 6.4 A (top panel) and the schematic 
(001) ribbon model (bottom panel). The inset is schematic 
edge states for the whole Brillouin zone. 


the fundamental gap (denoted by E g ) has a different 
behavior. When L < 6.2 A, E g is negative, manifest¬ 
ing a semi-metallic nature. As L increases, E g first in¬ 
creases to a local maximum (~80 meV) at L = 6.4 A, 
and then decreases to zero at L c , forming a Dirac cone 
at X point as shown in the inset of Fig.Qb). The above 
strain-dependent band structure analysis demonstrates 
that monolayer SnTe without strain has an inverted gap 
at X point with a negative value. According to the par¬ 
ity criteria [15], it can not be Z2 topological insulator 
as there are two A' points in the whole BZ. However, 
bulk SnTe has rocksalt crystal structure, and therefore 
its (001) monolayer thin films possess z —> —z mirror 
symmetry which allows us to define the corresponding 
mirror Chern number N m [55] . Following the similar anal¬ 
ysis and argument El, we can infer monolayer SnTe is a 
2D TCI with non-zero mirror Chern number N rn = 2. 

The physical manifestation of non-zero mirror Chern 
number is the existence of gapless helical edge sates which 
are protected by 2 —► —z mirror symmetry. By per¬ 
forming DFT calculations on the 20nm-wide ribbon of 
monolayer SnTe, we explicitly simulated the topological 
protected edge states. As shown in Fig. [2] for the ribbon 
with small lattice constant 6.4 A in the inverted regime 
(Fig. 0b)), there are two bands with linear dispersion 
crossing each other and connecting from the conduction 
bands to valence bands. Due to time-reversal symmetry, 
another pair of Dirac-type edge states is also present in 
the gap as shown in the inset of Fig. [2] which is con¬ 
sistent with the previous analysis N m = 2. Since the 
bands forming the Dirac point have different eigenval¬ 
ues of z —> —z mirror operation, both inter- and intra- 
back scattering are forbidden between those states. Thus, 
if the mirror symmetry is strictly preserved, those edge 
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states can realize quantized transport. Moreover, the 
mirror is connected with the spin in the system, and 
transport current should be also spin-polarized. 

We also performed DFT calculations for multi-layer 
thin films using the experimental lattice constant of 
bulk SnTe, L exp = 6.312 A, complementing our previ¬ 
ous work [TIT. Figure [3]' a) reveals a non-monotonically 
change in E* with increasing thickness. Two critical 
thickness values divide the phase diagram into three 
parts. More specifically, except three, five, and seven lay¬ 
ers of SnTe thin films, all other thin films have negative 
(or, inverted) E*. For thin films of more than 17 lay¬ 
ers, E* almost stays constant, reaching the gap of bulk 
SnTe. As the thickness gradually decreases, the confine¬ 
ment effect increases, driving E* to decrease and reach 
zero at about 7 layers. However, the gap didn’t continue 
its monotonic change for thin films of less than 5 layers. 
Instead, gap closes again between monolayer and 3 layers 
due to the aforementioned negative inverted gap E* in 
monolayer SnTe. The thickness-dependent inverted gap 
is consistent with our previous work m in the multi¬ 
layer case, providing a complete phase diagram shown in 

Fig- Ufa). 

To elucidate the underlying mechanism of band inver¬ 
sions, we analyzed the orbital evolution in detail and the 
result is presented in Fig. [3](b) . Previous works show that 
the fundamental gap of IV-VI semiconductor is around L 
(X) points for bulk (thin films), and the bands near the 
Fermi level are mainly composed of p orbitals of cation 
and anion (Sn and Te for SnTe) [0 IS . As shown in 
Fig. |3][b), we start from the atomic limit where the chem¬ 
ical bonding between different atoms are ignored and fo¬ 
cus on the p orbitals. After taking SOC into account, the 
six-fold degenerate p orbitals split into four-fold degen¬ 
erate J = 3/2 states and two-fold degenerate J = 1/2 
states. J = 3/2 states are pushed up while J = 1/2 
states are lowered down in energy as shown in Step (i). 
We then included the hopping between those atomic or¬ 
bitals (chemical bonds). For the bulk states, the pp hop¬ 
ping between cation and anion vanishes at L points, and 
the sp and pd hopping will play a vital role and lead to 
the four-fold degenerate J = 3/2 states to further split 
into J z = ±3/2 and J z = ±1/2 states (z along [111] di¬ 
rection). More importantly, the splitting is even strong 
enough to lower Sn’s |J = 1/2, J z = ±1/2) states than 
Te’s |J = 3/2, J- = ±1/2) states, driving SnTe into 3D 
TCI phase (Step (ii)). While for PbTe, the splitting in¬ 
duced by chemical bonding is not large enough to pro¬ 
mote the band inversion although SOC is even stronger. 

For the thin films, the states around X points mainly 
inherit from the bulk states at L point. When films are 
thick enough, all the properties of those states such as 
orbital weights and band orders are preserved, which 
ensure all the thin films thicker than a critical thick¬ 
ness (7 layers for SnTe) have negative inverted gap at 
X points and possess the nontrivial topological phase, 
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FIG. 3. (a)Thickness-dependent band gap at X point of SnTe 
thin films. Except thin films with three, five, and seven atomic 
layers, all other thin films have negative inverted gap and thus 
nontrivial topological phases \Nm\ = 2. (b) Schematic of 

band evolution in SnTe. (i) in the atomic limit, spin-orbital 
coupling splits the p orbital; (ii) in the bulk limit (z is set 
along [111] direction), strong chemical bonding between Sn 
and Te induces band inversion between | J = 3/2, J z = ± 1/2) 
of Sn and |J = 1/2, J z = ±1/2) of Te, which drives bulk 
SnTe into 3D nontrivial TCI phase; (iii) for SnTe ultrathin 
films of 3, 5, and 7 layers, the confinement effect surpasses 
the inverted bulk gap and removes the band inversion formed 
in Step (ii); (iv) for monolayer (z is set along [001] direction), 
very different chemical bonding environment of p x ,y and p z 
leads to the asymmetric crystal field effect. As a result, the 
onsite energy of p z will be lower (higher) than p XiV for Sn 
(Te), pushing all the \J Z = ±1/2) states from Sn (Te) to lower 
(higher) energy and hence re-introducing band inversion. 

i.e. \N m \ = 2. However, as the thickness decreases, the 
confinement effect will increase, hence enlarge the energy 
gap between anion and cation like in the normal semi¬ 
conductor, giving rise to vanishing band inversion (see 
Step (iii) in Fig. [3])b)) and topological phase transition 
as shown around 7 layers in Fig. [3])a). 

So far, all p x ,y,z orbitals are treated equally. This is 
reasonable for bulk and thick films, since the difference 
between all p x ,y,z orbitals is negligible. However, it is no 
longer true for atomically thin 2D materials as the chem¬ 
ical environment (crystal field) for p z and p x , y is com¬ 
pletely different, where 2 is along [001] direction. For 
Px,y orbitals, the crystal field is similar to the bulk where 
the p Xt y orbitals between anion and cation form both a 
and 7T type bonding. However, a bond is absent for the 
p z orbitals in monolayer. Such kind of asymmetric crys¬ 
tal field leads to the on-site energy of p z higher or lower 
than p XtV for cation (Sn) or anion (Te), respectively, i.e. 
Ap z (Sn) < 0 and Ap z (Te) > 0. This onsite energy differ- 
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FIG. 4. The band gap at X point (Ef) and the fundamental 
gap E g for all monolayer IV-VI semiconductors. All the ma¬ 
terials have negative inverted gap at X point and |Am| = 2. 
GeY (Y=S, Se, Te) are semi-metal and all other materials are 
fully gapped semiconductors. Among them, monolayer PbTe 
has the largest fundamental gap of 260 meV. 


ence Ap z will push down (up) all the J z = ±1/2 states 
composed of Sn cation (Te anion) orbitals since all those 
states couple with each other and involve with L z = 0 
states ( i.e. p z orbital). From the DFT calculations, we 
found the asymmetric crystal field effect is strong enough 
to invert the bands between anion and cation again as 
shown in Step (iv) in Fig. |3](b). We also confirmed this 
band inversion mechanism and the corresponding topo¬ 
logical phase transition by tuning the onsite energy dif¬ 
ference in TB calculations. 

It is worth to emphasize that the above band inversion 
mechanism in atomically thin films doesn’t depend on the 
band inversion in bulk. Instead, it critically relies on the 
competition between the trivial band gap and the strong 
crystal field effect Ap z on the p x ,y,z orbitals in the mono- 
layer. If Ap z is big enough, we can also achieve nontrivial 
topological phases in monolayers of other materials even 
if their bulk forms don’t have inverted gap. Motivated 
by this mechanism, we calculated other monolayer IV- 
VI semiconductors XY (X=Ge, Pb Sn and Y=S, Se, Te) 
and found that all them are 2D TCIs with |JV TO | = 2 
in the optimized structure. The detailed band struc¬ 
tures can be found in the Supplementary Materials. In 
Fig. 4j we plot the inverted gap at X point (£//) and 
the fundamental gap E g . Except GeY (Y=S, Se, Te) 
which are semi-metals[52| with the largest inverted gap 
at X , all other materials including PbY and SnY (Y=S, 
Se, Te) are fully gapped semiconductors. Among the 
latter ones, monolayer PbTe has the largest fundamen¬ 
tal gap of about 260 meV, comparable with all existing 
and proposed 2D topological materials. We have con¬ 
firmed nontrivial topological phases in all the materials 
by direct edge state calculations [13 33 . From the ma¬ 
terials’ perspective, PbTe is perhaps the most promising 


material to realize the novel application of topologically- 
protected edge states. First, PbTe holds the highest mo¬ 
bility among all those materials [47] . Second, PbTe can 
have both n-type and p-type carriers HI S3 and its qual¬ 
ity is easier to control in experiments m ■ In addition, 
except GeS, all other materials can realize the rocksalt 
structure under proper conditions [1HU5T| . Therefore, we 
expect that monolayer IV-VI semiconductors in rocksalt 
structure can be synthesized in experiments and realize 
the nontrivial topological phase presented above. 

CONCLUSION 

In conclusion, using first-principles calculations we the¬ 
oretically studied the electric structures of monolayer IV- 
VI semiconductors along [001] direction. Remarkably, 
contrary to the expectation from the quantum confine¬ 
ment effect, all these materials are 2D TCIs with mirror 
Clrern number |jV m | = 2, and monolayer PbTe possesses 
the largest fundamental band gap of 260 meV. The ID 
metallic edge states in these monolayer TCIs are pro¬ 
tected by top-to-bottom mirror symmetry and robust 
against any perturbation preserving the mirror symme¬ 
try. The unexpected nontrivial topological phase origi¬ 
nates from the strong asymmetric crystal held effect on 
p x<y and p z orbitals due to the unique characteristic of ul- 
trathin 2D materials. Our work makes monolayer IV-IV 
semiconductors a promising materials platform to realize 
novel applications of topological edge states, and offers 
a new strategy to find other 2D topologically nontrivial 
materials. 

METHODS 

The calculations were performed in the framework 
of first-principles density-functional theory (DFT) im¬ 
plemented in the Vienna ab initio simulation package 
(VASP) [53]. We used the generalized gradient approx¬ 
imation (GGA) of exchange-correlation function in the 
Perdew-Burke-Ernzerhof (PBE) form [53]. The projec¬ 
tor augmented wave method [55i was applied to model 
the core electrons. Monkhorst-Pack k-point sampling of 
10x10x1 was used for slab calculations. Energy cutoff 
of the planewave basis was fully tested for all the materi¬ 
als, and atomic structures were optimized with maximal 
residual forces smaller than 0.01 eV/A. Spin-orbit cou¬ 
pling (SOC) was included in all calculations. 
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(X=Ge, Sn, Pb and Y=S, Se, Te) in details. All those 
materials have negative inverted gap at X point, indi¬ 
cating that they are nontrivial 2D topological crystalline 
insulators with non-zero mirror Chern number \N m \ = 2. 
All the GeY (Y=S, Se, Te) are senri-nretals, while all 
other materials are semiconductors. The large fundamen¬ 
tal gap in PbY and SnY comes from the strong spin-orbit 
coupling of Pb and Sn. 

We define the sign of inverted gap E* in all monolayer 
IV-VI semiconductors XY using the strain effect. We 
have also confirmed the nontrivial topological phases in 
all these materials by directly calculating the edge states. 


SUPPLEMENTARY MATERIALS 


In the Supporting Information, we show the band 
structures of all monolayer IV-VI semiconductors XY 


E (eV) E (eV) £ (eV) 


7 



FIG. 5. Band structures of all monolayer IV-VI semiconductors XY (X=Ge, Sn, Pb and Y=S, Se, Te). 




































